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Double-polarization observables in the reaction ep —>■ e'j)''y have been measured at = 

0.33 (GeV/c)^. The experiment was performed at the spectrometer setup of the A1 Collabora¬ 
tion using the 855 MeV polarized electron beam provided by the Mainz Microtron (MAMI) and a 
recoil proton polarimeter. From the double-polarization observables the structure function is 
extracted for the first time, with the value ( — 15.4 ±3.3(stat.) -2 4 (syst.)) GeV“^, using the low-energy 
theorem for Virtual Compton Scattering. This structure function provides a hitherto unmeasured 
linear combination of the generalized polarizabilities of the proton. 

PACS numbers: 13.60.Fz, 14.20.Dh, 25.30.Rw 


INTRODUCTION 

Polarizabilities parametrize the response of systems 
composed of charged constituents to electric and mag¬ 
netic external fields. For the proton they contain in¬ 
formation about the QCD interaction in the very low 
momentum-transfer domain where the coupling constant 
ttstrong diverges. Since no static field of sufficient strength 
can be produced experimentally they are measured by 
means of Real Compton Scattering (RCS). Now due 
to the availability of powerful electron accelerators also 
Virtual Compton Scattering (VCS) can be investigated. 
VCS allows for the determination of Generalized Polariz¬ 
abilities (GPs) as function of the initial photon virtually 
as first pointed out in [ij for atomic nuclei and in [ij 
for nucleons. Just as the form factors Ge and Gm give 
access to the spatial density of charge and magnetization 
in the nucleon, the GPs give access to such densities for a 
nucleon deformed by an applied quasi-static electromag¬ 
netic field Out of the six lowest-order GPs of the 

proton, the electric and magnetic GPs have already been 
the subject of experimental investigation at MAMI d, Q , 
Bates @ and JLab [^. The four remaining ones, called 
the spin GPs, are still totally unknown experimentally. 
This Letter presents the first measurement of a double¬ 
polarization observable in VCS, with the aim of gaining 


insight into the spin-GP sector of the nucleon for the first 
time. 


FORMALISM AND NOTATION 

VCS is experimentally accessed through the photon 
electroproduction reaction ep —)• e'p'j. At low energy it 
can be decomposed into a dominant Bethe-Heitler (BH) 
part, a VCS Born (B) part and a VCS non-Born (nB) 
part, as shown in Fig.[TJ The contributions of the Bethe- 
Heitler and Born processes (BH-f B) can be exactly calcu¬ 
lated using as input only the form factors of the nucleon. 
The non-Born part is parametrized at the first order in 
the real photon momentum q' by six GPs. With an un¬ 
polarized cross section measurement only two linear com¬ 
binations of the GPs can be determined. For extracting 
all the GPs, double-polarization measurements are re¬ 
quired. In this experiment, the beam-recoil polarization 
asymmetries were measured in the reaction ep —>■ e'p'"f. 

The main kinematical variables are defined in the (yp) 
center-of-mass (CM): the modulus of the momentum of 
the virtual photon gem, of the outgoing photon and 
the polar angle 9^-^ between the two photons. The vir¬ 
tual photon polarization e and the angle p between the 
leptonic and reaction planes complete the kinematics. 
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FIG. 1: Contributions to the photon electroproduction ampli¬ 
tude. The VCS non-Born is parametrized by the CPs, while 
the BH-bVCS Born contribution (BH-(-B) contains no CP ef¬ 
fect and is entirely calculable in QED. 


The Low Energy Theorem (LET) for the double¬ 
polarization observables was developed in imi and is 
only briefly recalled here. Experimentally, the double¬ 
polarization observable is determined via 


d^a{h, i) — d^a{h, —i) 
d^aih, i) + d^a{h, —i) ’ 


( 1 ) 


where i = x, y, ^ is the CM axis for the recoil proton 
polarization component, h = the beam helicity and 
d^a{h, i) the doubly polarized {ep —^ e'p'y) cross section. 
The LET expansion, which is valid below pion threshold, 
leads to: 

2dPa ’ ^ ^ 

where Ad®(T®^+® is the difference of the doubly polarized 
cross sections i) — —z) and d^a 

is the unpolarized (ep e'p'j) cross section. (4>q(.-^) is 
a phase-space factor. The non-Born terms AAf'^® are 
linear combinations of the VCS structure functions 
Pxtj Pttj Pl/ti PlT’ Plt^ which can be expressed as 
linear combinations of the six GPs. In particular, P^t is 
a linear combination of the structure functions Pll and 
Ptt, where Pll is proportional to the electric GP, and 
Ptt is a combination of two spin GPs: p(A^i,Mi)i 
p(Li,M2)i^ the latter corresponding to ^eim2 in the RCS 
limit of 0. Eor more detailed formulas, we refer 

the reader to refs. B13- 


EXPERIMENT 


The experiment was performed at the spectrometer 
setup of the A1 collaboration at MAMI [llj and details of 
the analysis can be found in [13, [13. Table U summarizes 
the two kinematical setups of the experiment. 


TABLE I: Parameters of the spectrometer setups: p is the 
central momentum and 6 the in-plane angle. Both settings 
are centered on the nominal kinematics defined by gem = 
600MeV/c, e = 0.64, q(.^ = 90MeV/c and p = 180°. They 
differ in the covered in-plane CM-angle 9~fj. 


Setup 

Beam 

Spectrometer A 

Spectrometer B 


E 

(MeV) 

Pproton 

(MeV/c) 

^proton 

^electron 

(MeV/c) 

^electron 

VCSQOa 

855 

620 

34.1° 

546 

50.6° 

VCSQOb 

855 

645 

38.0° 

539 

50.6° 


The polarized electron beam was delivered by MAMI 
with an electron energy oi E = 854.6 MeV and a lon¬ 
gitudinal beam polarization of Pb = 70% in average. 
The beam polarization was determined by a Mpller- 
polarimeter and was flipped on a random basis with 1 Hz 
on average to avoid false asymmetries. A beam current 
of 22 p,A was directed on a liquid hydrogen target with a 
length of 5 cm. The beam was rastered across the target 
to avoid local boiling. 


Two particles were detected in coincidence: the scat¬ 
tered electron in spectrometer B with a solid angle of 
5.6 msr and a momentum resolution of 10““^, and the re¬ 
coil proton in spectrometer A with a solid angle of 21 msr 
and the same momentum resolution. Thanks to the good 
timing resolution, of 0.9 ns (EWHM) for the coincidence 
time, no further particle identification was necessary. Be¬ 
hind the focal plane of spectrometer A, a proton po- 
larimeter determined the transverse components of the 
proton polarization in the focal plane (see refs. 14, l3|l. 


The reaction was further identified by the missing mass 
squared, i.e. the squared mass Mj. of the missing parti¬ 
cle X in the (ep e'p'X) process. The distribution 
shows a clean peak at zero from photon electroproduc¬ 
tion, which is well separated from the pion peak from 7r° 
electroproduction. 


Eor the analysis, events with a missing mass squared 
of —lOOOMeV^/c^ < < 4000MeV^/c"'^ and a coinci¬ 
dence time of —1.5 ns < < 1.5 ns were accepted as 

VCS events. Events from the side bands of the coinci¬ 
dence time distribution were used to estimate the back¬ 
ground contribution due to random coincidences. A cut 
was also required to eliminate the target endcaps. 

Eor the determination of the polarization observables 
only events below the pion threshold were selected, by 
a cut of < 126MeV/c. A standard set of further 
cuts were applied to ensure a clean reconstruction within 
the acceptance of the recoil polarimeter and to select the 
region of large analyzing power of the polarimeter [l3| . 
A sample of about 77 000 VCS events survived the cuts. 
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BEAM-RECOIL POLARIZATION ANALYSIS 


With the polarimeter, for each event the direction of 
the secondary scattering process in the carbon analyzer 
was determined. This direction is given by the polar and 
azimuthal scattering angles 0s and $s- The distribution 
of events is given by: 

a{Qs,^s,Ep) = ao[l + hPb Ac{&s,Ep)Py'P cos 

-hPbAcies,Ep)P^P sin^s], (3) 


it depends on the known anal yzing p ower of the car¬ 
bon analyzer Ac{&s,Ep) (see 12], [l^) and the trans¬ 
verse components of the proton double polarization ob¬ 
servable in the focal plane, Py^ and P^^. For a given 
set of CM polarizations the focal plane transverse 

components P^^ and Py^ can be calculated by Lorentz 
transformation, rotation and ray-tracing of the spin pre¬ 
cession in the magnetic field of the spectrometer. Thus, 
the CM polarizations can be fitted to the distribution of 
the azimuthal angle $s by a standard maximum likeli¬ 
hood method. This is the first step of the analysis. 

In principle, the statistical ensemble contains the in¬ 
formation for all three CM components of the polar¬ 
ization, since events with different orientation of the 
scattering plane have different paths in the magnetic 
field of the spectrometer, resulting in different trans¬ 
verse components in the focal plane. A detailed sim¬ 
ulation showed, however, that the longitudinal compo¬ 
nent P™ cannot be reconstructed with sufficient resolu¬ 
tion. Therefore this component was fixed in the analysis 
for each event to the value given by the BH-f B calcula¬ 
tion, i.e. P™ = Ad5cr®H+®/2d5crBH-tB^ simulation 
showed that this choice was sufficient to provide a non- 
biased fit of P™ and P™. A more realistic choice, i.e. 
adding a CP effect in the constraint on P™, was consid¬ 
ered only to evaluate systematic errors. 

The maximum likelihood fit yields the CM polarization 
components P™ and P™. The ht is made separately in 
five 9.yj bins to have sufficient statistical significance per 
bin. 


Tabled] summarizes the results. The obtained values 
for Py“ are compatible with zero within the uncertain¬ 
ties; this is consistent with the requirement that P™ has 
to vanish in strict in-plane kinematics. Globally, P™ 
has a negligible sensitivity to the GPs and almost all 
the new information is carried by P™, through the term 
which is of the form: 


A7W"®(h,x) = h {a^Ptr + a^ExT 

+o.^Ptt + ) ( 4 ) 

with af being known kinematical coefficients [13- 
Figure m displays the measured P™ component as five 
solid points. These points have been projected to the 


TABLE II: Results for the double-polarization observables. 
Pj:™(raw) and Pjr™(proj.) are the fitted Pjr™ component be¬ 
fore and after the projection to the nominal kinematics, re¬ 
spectively. APjr™ (stat.) is the statistical error on Pj:™(proj.), 
while AP™ (syst.) are systematic errors (see text). Negative 
6j.y values are conventional for ip = 180°. 



9.y.y 

-170° 

-150° 

-130° 

-110° 

-90° 

pcm 

(raw) 

0.047 

0.012 

-0.043 

0.020 

-0.020 

AP™ 

(stat.) 

±0.066 

±0.053 

±0.038 

±0.041 

±0.050 

pcm 

(raw) 

-0.220 

-0.269 

-0.215 

-0.177 

-0.067 

pr 

(proj.) 

-0.209 

-0.257 

-0.201 

-0.142 

-0.041 

AP°“ 

(stat.) 

±0.049 

±0.040 

±0.030 

±0.027 

±0.027 

AP°"' 

(syst.l) 

±0.001 

±0.011 

±0.007 

±0.009 

±0.004 

AP°“ 

(syst.2) 

±0.030 

±0.001 

±0.003 

±0.020 

±0.030 

AP°” 

(syst. 3) 

±0.010 

±0.020 

±0.020 

±0.020 

±0.010 


nominal kinematics (gcm,e, ^cmi ^-^d of Table |T|, com¬ 
pleted by the values of 0-^.^ of Table (Til This projection is 
based on the expected LET behavior of the polarization 
observables as a function of the kinematics. 

The statistical error is provided by the fit. Systematic 
errors on P°™ have been determined as coming from: 1) 
a beam polarization uncertainty of ± 1.2% ; 2) changing 
the constraint on P™; 3) uncertainties in the kinematical 
projection. Other systematic effects, due to instrumental 
asymmetries in the proton polarimeter or due to random 
coincidences under the time peak were found to be neg¬ 
ligible. 


STRUCTURE FUNCTION ANALYSIS 

As a next step, a ht was performed with the aim of 
determining individual GPs (including the spin GPs). 
The principle is again to use the likelihood method, this 
time fully unbinned. The non-Born terms AAi’^^{h,x) 
and AA4’^^(h, y) in the numerator of P™ and P™' were 
replaced by their analytical expression in terms of the 
GPs M- The cross section d^a in the denominator of 
Pjr™ and P™ was hxed to its value given by the unpo¬ 
larized LET expression, using our previously measured 
structure functions {Pll — Ptt/^) and Plt [3]- As an 
outcome, it turned out that the data were not precise 
enough to extract individual GPs. However, if one uses 
structure functions, i.e. combinations of GPs, instead of 
GPs directly, one gets a signihcant result for P^j. as we 
show in the following hnal step of the analysis. 

The unbinned maximum likelihood method is again 
used. The non-Born terms AAI"®(/i, x) and AAl”®(h, y) 
are replaced by their analytical expressions in terms of 
the structure functions, as e.g. in Eq. O- The denomi¬ 
nators of Pi:™ and P.'r™ are treated as above. 

oL y 

Exploratory fits showed that AAI"® is sensitive mainly 
to PxT^ among the four structure functions entering 
Eq. (|3]). Therefore, the other three: P^y,Pyy and Pjjf, 
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TABLE III: The complete set of VCS structure functions at = 0.33(GeV/c)^ as calculated by two models: (I)= DR 
formalism [l^ (with Aa = 1.80 GeV, Ap = 0.75 GeV and the MAID03 version), (II) = HBGhPT at 0{p^) [l3], including the 
7r°-pole term (or anomaly). Only six of these nine structure functions are independent, e.g. the ones in the first six columns. 
Pttj Pttj and Pijr are fixed to the values of this Table when fitting (see text). All calculations are done with the proton 
form factors of ref. [3 • This also holds for Table [V] 

Model Structure Functions (GeV“^) 

model Ptt Ptt Plt Plt Plt Pll Ptt Plt Put 

(I) DR model 0.97 -0.44 -10.83 -1.43 -2.43 22.40 -1.58 -1.34 -1.21 

(II) HBGhPT 0{p^) 2.05 0.62 -10.57 -4.21 -5.34 15.07 -6.89 -3.03 -0.86 


cannot be fitted. However, their influence can be inves¬ 
tigated by inserting several model predictions and fitting 
only. This implies a model dependence of the ex¬ 
tracted results, but we show in the following that it is 
relatively small compared to the statistical uncertainty. 

The structure functions that need to be fixed are 
only Pr^rp^Pj,^ and i.e. the ones appearing in 

AM''^{h,x) and except P^x- The max¬ 

imum likelihood fit was done with three rather different 
assumptions for these hxed structure functions. In fit 
“I” they were set to values calculated by the Dispersion 
Relation (DR) model [l^ , cf. the first line of Table IIIII 
In fit “H” they were set to values calculated by Heavy 
Baryon Chiral Perturbation Theory (HBChPT) [la, cf. 
the second line of Table Hill In fit “HI” they were all set to 
zero. These different choices lead to the following results: 

= -15.4,-17.7 and -14.1 GeV"^ for fits “F, “H” 
and “HI” respectively. We consider fit “I” as the central 
one, yielding our final result for Pp^i the two other 
results are used to estimate the model-dependent error. 

The statistical error on Pp^ is provided by the max¬ 
imum likelihood fit. The systematic error comes from 
several main sources, which are estimated in Table lYl 
The first contribution is obtained by changing the beam 
polarization by ± 1.2% in the analysis. The second con¬ 
tribution is estimated by performing the fit with several 
form factor parametrizations [l9l - l^ ; the maximal spread 
of the results gives the magnitude of the error, which re¬ 
mains small. The third contribution is related to the 
treatment of the error is obtained as the difference 
in the fitted result when we fix P™ to its BH-I-B value, 
or when a GP effect is added to it. The fourth contribu¬ 
tion is due to model dependence; it is determined from 
the differences between the various fits (“H” — ”1” and 
“HI” — ” I”). In Table IIVI each partial systematic error 
has been symmetrized except the fourth one which is the 
largest and most asymmetric. The total systematic er¬ 
ror is calculated as the quadratic sum of the errors of 
Table US for each sign separately. 

Our final result for Pp^ is presented in Table 0 It is 
compared to theoretical values from HBGhPT and DR 
calculations. The absolute value of the result is larger 


TABLE IV: Systematic errors in the extraction of Ppx 


Error Type 

Error Value (GeV-"") 

beam polarization (±1.2%) 

± 0.53 

proton form factors 

± 0.10 

constraint on P™ 

± 0.47 

fixed structure functions 

±1.26/ - 2.29 

Total systematic error 

±1.45/ - 2.40 


than in most theoretical calculations. Some features of 
the models are worth noting: In HBGhPT some of the 
GPs have a bad con verg ence with respect to the order 
of the calculation ^ 241, and this may affect the model 
value of P^T- DR model the spin GPs are entirely 

fixed, but the scalar GPs contain an unconstrained part 
that has to be fitted from experiment. In particular P^^ 
depends, via the structure function Ppp, on the free pa¬ 
rameter Aq, which determines the electric GP. Table El 
shows this dependence for a realistic range of values for 
Aq,. We note that the DR model has a lower limit for 
Ppx of —13.1 GeV“^ (for A^ = oo). 


A graphical representation of our result is shown in 
Fig. [2l The central solid curve is obtained by calculating 
the polarization component P™ at the nominal kinemat¬ 
ics, based on Eqs. @ and The calculation uses the 
results of fit “I” (see above), i.e. Pp^ = —15.4GeV“^ 
and the other three structure functions set to their DR 
value of Table Hill Using the results of fit “H” instead of 
“I” yields a very similar curve. The deviation from the 
BH-I-B calculation (dashed curve) is a clear signature of 
the polarizability effect. 


In conclusion, we have measured for the first time 
double-polarization observables in VCS from the proton 
below the pion threshold. The analysis was based on the 
theoretical formulation of the LET for polarized VCS, 
and the experimental use of recoil proton polarimetry. A 
clear polarizability effect was observed in the P“™ po¬ 
larization component. We extracted one new structure 
function, P^tj ^^.d found a value that is larger in magni¬ 
tude than most theoretical calculations. Therefore, this 
measurement provides a valuable and entirely new con- 
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TABLE V: Our measured value of and several model 
predictions at = 0.33 (GeV/c)^. For the DR model, the 
(a), (b), (c) cases correspond to different values of the Aa 
parameter: 0.6, 1.2, and 1.8 GeV respectively (see text). 


PtT (GeV-") 

This experiment 

15.4 zb 3.3(stat.) —2.4 (syst.) 

DR model 

-3.7 (a) , -8.7 (b) , -10.8 (c) 

HBChPT 0{p^) [17] 

-10.6 



FIG. 2: Measured recoil proton polarization component 
in the CM frame. The five points with their statistical error 
are the result of the first-step fit. The solid curve is calculated 
using our result forP^j. (see text); the shaded band represents 
the statistical uncertainty. The dashed curve is the BH-fB 
calculation of P)r™, i.e. without any GP effect. 

strain! for models of nucleon structure; although it does 
not allow one to further disentangle the scalar and spin 
GPs of the proton. 
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